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^ ' Abstract 

O ' 

Universality in unitary invariant random matrix ensembles with complex matrix 
elements is considered. We treat two general ensembles which have a determinant 
factor in the weight. These ensembles are relevant, e.g., for spectra of the Dirac 
operator in QCD. In addition to the well established universality with respect to 
the choice of potential, we prove that microscopic spectral correlators are unaffected 
when the matrix in the determinant is replaced by an expansion in powers of the 
matrix. We refer to this invariance as logarithmic universality. The result is used in 
proving that a simple random matrix model with Ginsparg- Wilson symmetry has 
the same microscopic spectral correlators as chiral random matrix theory. 
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1 Introduction 



Random matrix theory (RMT) deals with universal spectral properties of en- 
semble averages of large random matrices p] . This universality is revealed in 
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the microscopic limit. In the microscopic limit we consider the spectral prop- 
erties of eigenvalues A of matrices of dimension N in the limit as iV — > oo 
with x = NX fixed 0. Studies of universality in RMT can be divided into two 
categories. In one category the models studied have unitary invariance (see, 
e.g., |J) while in the other this invariance is violated (see, e.g., [[|). In this 
paper we consider the former category. The partition function defining the 
ensemble, hence, can be expresses in terms of the eigenvalues of the random 
matrices rather than in terms of the matrices. We consider random matrix 
models where the partition functions involves a general potential V^(A). 
This general potential is a sum of a regular and a logarithmically singular 
part. This paper introduces and proves logarithmic universality. Logarithmic 
universality expresses the invariance of microscopic spectral correlators with 
respect to deformations of the logarithmically singular part of V^ Q )(A). The 
logaritmically singular part in the eigenvalue representation originates in the 
determinant of a matrix D. 

Two random matrix ensembles will be studied: The unitary ensemble (UE) 
relevant, e.g., for QCD in 3 dimensions |5|-|7] and the chiral unitary ensemble 
(%UE) relevant, e.g., for QCD in 4 dimensions @,§,f|. In the applications of 
RMT to QCD the matrix, D, in the determinant is analogous to the Dirac 
operator. For the xUE the logaritmic universality includes deformations which 
violate both the chiral symmetry condition {D, 75} = and the hermiticity. As 
an example we will show that a simple random matrix model which satisfies 
the Ginsparg- Wilson relation {_D,7 5 } = -D75-D ]T0| shares the microscopic 
correlators of the xGUE. 

In the first part of this paper, we introduce and prove logarithmic universal- 
ity in the unitary ensemble. In order to do this we make use of a method 
established by Kanzieper and Freilikher [II]. In the second part, we introduce 



logarithmic universality in the chiral unitary ensemble and use the universal- 
ity of the UE to prove the logarithmic universality of the microscopic spectral 
correlators of the %UE. Finally, we consider the example with the Ginsparg- 
Wilson symmetric model. 



2 Logarithmic universality in the unitary ensemble 

The unitary ensemble is defined by the partition function (see e.g. fl 



Z VE = J dMdet 2a Me- NTrV M (1) 
00 N 

~ / II (dA i A l 2a ex P {-W(A i )}) |A(A)| 2 , (2) 



-00 i=1 
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where the Aj are the eigenvalues of the hermitian N x N matrices M and 
V(M) is an expantion in even powers of M 

p 

V(M) = J2^M 2k , £ 2fc G M with £ 2p > . (3) 

k=l 



The Vandermonde determinant, 

N 

A(A) = J] i\ " A,) , (4) 

i>i=i 



is the Jacobian describing the change to angular coordinates Finally, the 
measure of the matrix integration is the Haar measure. 

The hard edge, bulk, and soft edge correlators of this ensemble have been 
studied in detail and have been shown || to be independent of the coefficients 
g2k in the potential V in the microscopic limit. 

We now introduce a new set of deformations and show that the microscopic 
spectral correlators are invariant with respect to these deformations. 

We shall prove that the microscopic spectral correlators of the partition func- 
tion 

Z VE = J dMdet 2 " (/(M)) exp{-NTtV(M)} (5) 

oo N 

- I n(dA,:[/(A,)] 2Q exp{-W(A J )})|A(A)| 2 (6) 

-oo i=1 

are independent of the coefficients G M in the polynomial 

oo 

/(M) = ]T/ 2fc+1 M 2fc+1 (7) 

k=0 

with /i ^ 0. The expansion defining / is assumed to have suitable convergence 
properties. We denote this independence as logarithmic universality. 

In the proof of this assertion, we shall follow, and extend where necessary, 
the argument of Kanzieper and Freilikher [IT]. At present [11] provides the 



most far-reaching proof of the range of the universality class. At the same 
time, this method enables us to explore the spectral properties of the matrix 
model defined in eqs. @ and (|6|) for an arbitrary spectral range. Without loss 
of generality, we can choose fi = 1. 



3 



2. 1 Proof of logarithmic universality in the UE 



For notational simplicity, we will absorb the /-dependent part of the integrand 
in the partition function of eqn. (^) into a general potential 

VW(A) = V(A) + kL1(A) = t ff^ - ^log[/(A)] 2 • (8) 

k=l 



In order to determine the kernel, Kn(\, A'), from which all spectral correlators 
follow, 

p(Ai, A 2 , . . . , A fl ) = det K N (X a , A 6 ) , (9) 

Ka,b<s 



it is sufficient to find the set of polynomials P n (X) orthonormal with respect 
to the measure d/i(A) = dAu;(A) = dAexp{— NV^ a '(X)} on the real axis 



d//(A)P„(A)P m (A) = 8 mn . (10) 



Along with the polynomials P n (A), we will require the functions y n (A) 

^ n (A)EEP„(A)exp{-^V( a )(A)}, (11) 



which are by construction orthonormal on the real axis with respect to the 
measure dA. The kernel can be expressed in terms of the (p^ as 



JV-l 

k n (\,\>)=y,M^MX) (12) 

k=0 

Vjv(A')Vjv-i(A) - p N (\)<p N -i(\') 
= c N — x , (13) 



where the latter equality is derived using the Christoffel-Darboux formula fl4" 
The three term recursion relation 

AP n _i(A) = c n P n (A) - c n _iP n _ 2 (A) (14) 



determines the coefficients c n of which we have just encountered c^. The aim 
is to derive a differential equation for the orthonormal functions <^fc(A) ||11||. 
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In order to do so, we write the derivative of P n (A) as 
dP n (A) 



dA 



A n (A)P n _ x (A)-P n (A)P n (A) 



(15) 



where the functions A n (X) and B n (X) are |TT 



A n (\)=Nc n 



d/z(t) /dy^(t) dV^A)' 



-OO 
OO 



B n (X)=Nc n J 



t-X 



dt 



dA 



P n (t)P n ^{t) . (17) 



These expressions for A n (X) and B n (X) can be obtained by expressing the left 



hand side of eqn. ([Tq) through an expansion 



dP n (A) 7 m dP n (t)r^ 

d M*) ^ 2^ Pfc(A)P fc (t) 



dA 



df 



(18) 



fc=0 



and then rewriting this integral using partial integration and the Christoffel- 
Darboux formula. 



Using eqs. ([[6]), (|T7|), and the three term recursion relation eqn. (fH|), one finds 
the following identity 



B n {X) + P n _x(A) + iV d ^ (A) = — A>-i(A) . 



dA 



Cn-l 



(19) 



We are now able to derive a differential equation for the orthonormal functions 



ip n . Differentiating eqn. ([Tq ) with respect to A and making use of the recursion 
relation eqn. (fH] ) and of the defining eqn. ([TT]) for (p n , we obtain \T1 



dWA) -^)^ + »„(A)=0 



dA 2 



dA 



(20) 



where 



T n (X) = — \ogA n (X) 
dA 



(21) 



and 
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g n (X) = — —A n (X)A n ^ 1 (X) - B n (X) + — — — (22) 
Cn-i \ 2 dA J 

d (n m + N dv(a) W \ dlog A n (X) ( N dV^(X) \ 

It is straightforward to obtain the kernel, provided that we can solve these 
equations for cp^ and <pn-i, cf. eqn. (p~3|) . 



In accordance with |TT[ we choose the following convenient representation of 
the functions A n (X) and B n (X), 



AnW=^l(X)+aA^ g (X) (23) 
B n (X) = BW(A) + aB<&(A) , (24) 

where the regular parts depend only on the polynomial terms in the general 
potential, eqn. (|8|), and the singular part depends on the term log[/(A)] 2 . Since 
we have not changed the polynomial part, we are led to the same expressions 
for the regular parts of A^(X) and -Bjv(A) as obtained in |TT|. As our interest 



is the microscopic limit, we evaluate the differential equation, eqn. (£(]), in 
this limit. In doing so it is essential to keep track of the iV-dependence of 
the various terms. We will assume that the recursion coefficients c n and the 
functions A^ and approach smooth functions in the large- N limit so 
that 



c N±1 = c N + 0(l/N) (25) 
A^=A^ + 0(l/N) (26) 
B^=B^ + 0(1/N). (27) 

With these conditions, the function A^ is related to the macroscopic spectral 
density [0,0 



p(X) = j™ 0O ^5 ) (A) v / l-(A/a) 2 , (28) 



where a = lmi/v^oo 2cn can be identified with the endpoints of the spectrum. 
From this, we infer that A^(X) is of order N in the limit iV — > oo, provided 
that A is sufficiently far from the spectral endponts. 



Using the assumption of eqn. ( [27D to evaluate eqn. (0) in the large- N limit, 
we obtain 



Bn{X) + — - 2^ A n{A) ~ 2 losing - ^sing) ■ < - > ) 
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The singular parts Ag in g(A) an d BlingW depend on /(A). We will now prove 
that the /-dependent terms are are smaller than the /-independent terms by 
at least a factor of 1/N in the microscopic limit. 

First, we determine the A dependence of A^(A) and -B s ^g(A). The singular 
part A^JX) reads 



A (N)(^_ „ f Mt) ( dV sin S (t) dK ing (A) \ 2 
4ing(A) - ~C N J j-^l — — | P N (t) 



-oo 
oo 



d^t) / Mlog[/(t)] 2 dlog[/(A)H 2 

~ CN 1 t^x[ — a* dA — r N{t) 



-oo 
oo 



- - 2CN J — { JWW J ■ (30) 

— oo 

We now introduce the defining expansion of /(A) and expand l/(i — A) in A: 



oo 



(/'(«)(* + AA 3 + ...)- /(t)(l + 3/ 3 A 2 + ...)) P 2 N (t) 



A + / 3 A 3 + . . . J " w ( t t 2 t /(f) 

— OO 



/ d ^ h + - t m + - \ p2 ^ (31) 



1 + / 3 a 2 + . . . j ^ w t t 2 t /(*) 

—00 

In the final equality we have used the fact that both the weight u(t) of the 
measure and P%{f) are even in t. P^(t) is even in t since P n (t) has parity 
Pn{—t) = (— l) n P n (t), as can be shown using the three term recursion relation 
eqn. (|T4]) iteratively. 



The singular part of -Bat (A) is evaluated in an almost analogous fashion. The 
only extra ingredient is that one must use the recursion relation eqn. ([14]) in 
order to perform the integrations. We find that the term proportional to 1/A 
is independent of /. In addition to this term, which is identical to the one 
obtained in [0]] , we obtain a constant term and terms proportional to positive 
powers of A. Denoting the sum of the additional terms by Tb{N, A), we have 



1 — ( — 1) N 

BZl(X) = + T B (N, A) . (32) 
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To complete the proof, we determine the differential equation in the micro- 
scopic limit. Inserting B^j,(\) into eqn. (p9|), we find for the large- N limit of 
eqn.flTJ) 

N dV^(X) A a 
~2 dA = 2cat ~~ A" 2 



^(A) + ~ T\ = 77— ^4jv(A) - (-1)"- - -(T B (N - 1, A) - T B (iV, A)) .(33) 



Inserting this into the defining eqn. ([23|) for C/ n , we can determine its asymp- 
totic behaviour for large N 



^(A) = ^(A)(l--) + + (-1) -_ (34) 

+ (-l)"a) - ^(T B (iV - 1, A) - Tb(N, A)) 2 

+f (T B (iV - 1, A) - Tb{N, A)) ( A^( A ) _ 2( _i)iv« + 
2 ^cjy A A N (\)J 

-~{T B {N-1,X)-T B {N,X)). 



In the microscopic limit, two terms in C?at(A) are dominant (i.e., of order A^ 2 ). 
To see this we consider the N dependence in the microscopic limit of the 
various terms involved. As noted, eqn. (|28|) implies that 



45(A) oc N. (35) 



From the evaluation of the A-dependence of v4g in g(A) and -B sing (A) eqn. fl3~T|) 
and eqn. (B3), we have 



Affi(A)oc0(l) (36) 
Bffl(A)oc(l- (-1)> + 0(1) • (37) 

Finally, the macroscopic density and, hence, through eqn. (|28| ) is inde- 
pendent of A for A ~ 0. We can thus use eqn. ([H]) to obtain 

^(A) oc 0(1) . (38) 
Keeping only terms of order A^ 2 in eqn. (|20"D results in 



dV„(x) f (A ( g(0)) 2 , (-l^a-a 2 ! , . n 



(39) 
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Here we have changed variables to x = XN according to the microscopic scaling 
limit and divided by iV 2 . The important feature is that there are no references 
to the expansion coefficients fk in eqn. ([/]). Making use of eqn. (pT3[), the kernel 



of the UE follows cf. ||11|| . Since every orthonormal functions which enters in 
this kernel is universal with respect to deformations in both the potential V 
and logarithm / we infer that the microscopic correlators obtained from this 
kernel by eqn. (P) are universal. This completes the proof. 



The restriction to an odd polynomial in eqn. (|7[) was necessary because the 
weight in the measure need to be even and the lowest order term in eqn. (|7|) 
has to be linear. Following the lines of the proof given above, it is simple to 
prove that: If one replaces f(Xi) in Z\je with /(|Aj|) then without affecting the 
microscopic spectral correlators one can introduce even powers in the definition 
of /. That is, microscopic spectral correlators of the partition function 

oo N 

Z - I II (dM/(|A,|)] 2a exp{-iW(A 4 )}) |A(A)| 2 . (40) 

-oo i=1 

are independent of the coefficients fk € K of the polynomial 

oo 

/(A) = E/^ (41) 

fc=i 



with /i ^ 0. 

While this partition function to our knowledge does not have any explicit 
matrix analogue it will nevertheless prove useful as we turn to the chiral 
unitary ensemble. 



3 Logarithmic universality in the chiral unitary ensemble 



Chiral unitary random matrix ensembles were originally introduced to 
model the consequences of chiral symmetry in the spectrum of the Euclidean 
Dirac operator of QCD. The structure of the random matrix representation, 
M, of the Dirac operator was determined by the vanishing anticommutator 
of M and 75. This anticommutation has one particularly simple spectral con- 
sequence: The eigenvalues come in pairs of opposite sign. In extending the 
universality of the chiral ensemble, we shall again replace the matrix, M, in 
the determinant by an expansion in powers of the matrix. However, we must 
be somewhat careful because we will need to retain a remnant of chiral sym- 
metry. Here we show that the condition for universality of the microscopic 
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spectral correlators is complex conjugation symmetry of the spectrum and 
not chiral symmetry as established through the anticommutation with 75. 

To this end, we we will prowe the assertion that: The microscopic correlators 
of the partition function 

Z xVE = J dMdet" (/(M)) exp{-NTrV(WW^)} , (42) 

where 6t is an integer, are independent of the coefficients fk G R of the poly- 
nomial 

00 

/(M) = £ hM k (43) 

k=l 

when /1 / 0. 

By definition the hermitian 2N x 2N matrix M is block off-diagonal: 

( iw\ 

M=\ , (44) 

\iWi J 

where If is a complex N x N matrix. The non-singular potential is given by 

v ~ 

V(WW*) = ^(WW^) k , g k G R with g p > . (45) 
fc=i ^ 

Again, it is sufficient to concider /1 = 1. 

3.1 Proof of logarithmic universality in the xUE 

The strategy of the proof is first to write the %UE partition function in an 
eigenvalue representation and then to relate this partition function to that of 
the eqn. (TO) introducing appropriate shifts of the arguments. This will allow 
us to express the kernel of the xUE in terms of the orthonormal functions of 
the UE and, hence, carry over the logarithmic universality of the UE to the 
XUE. 

We rewrite the M integration in the generating functional in angular coordi- 
nates W = VAW, where V and U are unitary matrices and A is a diagonal 



10 



matrix with real and positive entries Aj, i = 1, . . . , N. This is useful since 
-Z x ue is invariant under the transformation of W — > WWV. This invariance 
is obvious because the transformation of M and hence of f(M) is unitary. 



Under the unitary transformation of M, the measure dM changes to [13[] 

N 

dM = Y[d\ 2 i [dV}[dU}\A(\ 2 )\ 2 , where A(A 2 ) = J] (A, 2 - A 2 ) . (46) 

j=l i<jf 



The eigenvalues A 2 of WW^ 1 are related to the eigenvalues of M, which are 
±i\. 

It follows from the reflection symmetry of the spectrum of M and the choice 
of real coefficients fk that the spectrum of f(M) has complex conjugation 
symmetry The determinant factor can therefore be written 

N 

det(f(M)) = l[f(i\ j )f(i\ j y. (47) 

5=1 



Eqn. ( |4"7| ) show that the determinant is a function of the A 2 . Thus, the complex 
conjugation symmetry of the spectrum of f(M) is the remnant of chiral sym- 
metry which we have retained to ensure that the determinant can be written 
in terms of the eigenvalues of WW^ . 

Due to the cyclic property of the determinant and the trace, the integrations 
over V and U become trivial, producing irrelevant factors in -Z^ue- The par- 
tition function now reads 

oo N 

?xVE~ [ II (dA 2 |/(.A J )| 2fi exp{-W(A 2 )}) |A(A 2 )| 2 (48) 
o i= l 

OO 

~ / II (dz^llR^^M-NViz 2 )}) |A(z 2 )| 2 . (49) 
-oo i =1 



This partition function is not of the same form as that of the unitary en- 
semble eqn. (||). However, inspired by |12 |, we can carry over the universality 
argument of the preceding section. Given the set of coefficients fk defining the 
polynomial /, we introduce polynomials Pi^ a \z 2 ) orthonormal with respect 
to the measure dz\z\\f {iz)\ 2a ex.]i{—NV {z 2 )} on the real axis. Notice that it 
is possible to replace |A(z 2 )| 2 by a determinant involving only the even poly- 
nomials Pi^ a \z 2 ). The idea is to identify the polynomials P^' a \z 2 ) with the 
polynomials P%[ ,a \z) for the partition function eqn. (ffij). We will now show 
that this is justified ifa = a + l/2 and the expansion coefficients of / in 
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eqn. (|4l|) are related to those of / in eqn. (|43|) . This can be done provided that 
we identify 

V(z 2 ) = 2V(z) . (50) 



The proof is then completed by expressing the kernel of the xUE in terms of 
the orthonormal functions (fi of the partition function eqn. (fffip . Since these 
ortho normal functions are universal in the microscopic limit, the kernel is also 
universal in the microscopic limit. 

In order to find the set of coefficients {fk} in eqn. (f4lD such that |z||/(^)| 2a = 
jq z |-j2(q+i/2) - g sa ti s f] ec l f or a given set {/}, let us first expand the left hand 
side of this identity. Using the fact that fi = 1, we find 



Ml/»l 



2n 



E/ 2fc+ i(-i)^ 2fc+1 



+ (E(-i) fc+1 w^ +2 

\fc=0 



\k=0 



E\ f f I i \k+k' 2(k+k')+2 , f f I \fc+fc'+2 2(fc+fe')+4 
|_j2fc+l/2fc'+H-lJ ^ + J2fc+2j2fc'+2( 1 -lj ^ 



k,k'=0 



= £ (-l)^ ki+k * II [/2*,+i/»$+i + M+2/: 
= \ z \ 2 ^+h) + ... 



2fc^+2^ 



l 2 (Ef=i^+ fe D+ 2 «+ 1 



(51) 



Since only odd powers of |z| appear, it is appropriate to compare this expansion 
with that of /(M) 2( " +1/2) of eqn. (§!]) where f 2k +2 = 0. 



/ oo \ 2(fi+|) 

Lf(W)P a+ i> = (E/*fiN afc+1 



\fc=0 



fc ,fci,A:;,...,fc a ,K=0 j'=l 

:\z\ 2 ^ +1 ^ + ■■■ . (52) 



The proof that it is possible to choose the f 2 k+i for an arbitrary choice of 
coefficients fk such that the expansions of eqn. fl5"2]) and eqn. (|5T|) are identical 
to all orders can be made iteratively: First, note that the lowest-order terms 
in the two expansions are identical. Assume that the identity holds for each 
term up to power 2n + 2a + l of \z\. Now, the next term in the xUE expansion 
yields some coefficient for | 2 ;| 2 ( n + 1 )+ 2Q; + :1 - ) anc l this coefficient has to be matched 



12 



by 



k +k 1 +k' 1 + ...+k & +kl=n+l, fc,eN 3=1 



Since j\ = 1, we find that 2d + 1 of the terms in this summation are /2(n+i)+i- 
As /2( n +i)+i does not appear in any lower order term, one can choose it to 
and match the %UE expansion of eqn. (51) and the UE expansion of eqn. 



to order 2n + 2a + 1 in \z\. This completes the proof by induction. 

In order to determine the kernel of the xUE, we introduce orthonormal func- 
tions ipi as in eqn. fllT] ) 

#% 2 ) ee \z^\f(iz)\ & eM-jV(z 2 )}P[ f ' & \z 2 ). (53) 



Using the relation between the potentials imposed in eqn. (]50f) , we see that 

tf M & (54) 

The kernel is obtained by evaluating 

N-l 

i^W) = E #% 2 )#*V) • (55) 



Following the derivation in section 3 of [jX2[ , i.e., using the Christoffel-Draboux 
formula, we find 



(/(/),«+!)/ \ (/(/),«+§)/ \ (/(/),s+§), N (/(/>+§) 

K N (z , W ) = C 2 N 



z 2 — w 2 



It is immediately seen that the kernel is universal in the microscopic limit 
since the orthonormal functions tpi are independent of both V and / in that 
limit. 
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4 Example: The Ginsparg- Wilson relation 



As an indication of the utility of the results obtained above, we consider the 
partition function eqn. (f42|) with the choice 

1 00 M k 

D = f(M) = -[1 - e aM ] = - £ — a k ~ l , (57) 
a k=1 ft! 

where M is given by eqn. ([44]) and a is a real number. This choice falls within 
the scope of the results in section 3 since the expansion coefficients are real. 
Hence, the microscopic spectral correlators in that case are identical to those 
obtained for f{M) = M which is used when considering QCD in 4 dimensions 
with xUE 0,^3 . That model contains SU(a) x SU(a) chiral symmetry since 

{M, 75 } = 0, (58) 



where 

75 = diag(l,l,...,l,-l,-l,...,-l) . (59) 

It is readily seen that = 75-D75 and that the eigenvalues of D lie on a circle 
in the complex plane. Further, using eqn. ([58]) it is easy to show that 

{D, l5 } = aD l5 D . (60) 



This anti-commutation relation is known as the Ginsparg- Wilson relation [K| 
and implies a 577(a) x 577(a) symmetry of the model flT5| . This example shows 
that genuine chiral symmetry and the Ginsparg- Wilson-Liischer symmetry 
lead to identical microscopic correlators. The Ginsparg- Wilson relation is of 
relevance when constructing lattice actions without fermion doubling, see, e.g., 



[ 16 1 . Regarding a as the lattice spacing, it is clear that the Ginsparg- Wilson- 
Liischer symmetry will reduce to genuine chiral symmetry in the limit a — > 0. 
The present argument show however that the microscopic spectral correlators 
will be obtained for all a and that, in this regard, there is nothing special 
about the a — > limit. The first study of random matrix correlations in a 
Ginsparg- Wilson symmetric lattice action have appeared recently 



5 Summary 



In this paper we have extended the universality class of the unitary invari- 
ant random matrix ensemble and of the chiral unitarily invariant random 
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matrix ensemble to include logarithmic universality. Logarithmic universality 
expresses independence with respect to deformations of the logarithmic singu- 
lar part of the general potential. In the case of the xUE, this means that the 
original matrix structure in the determinant can be deformed by means of a 
series expansion in the original block matrix without altering the microscopic 
spectral behaviour provided that the spectrum associated with this series has 
a complex conjugation symmetry. This result was used to demonstrate that 
the matrix appearing in the determinant of the %UE partition function can 
be chosen to satisfy the Ginsparg- Wilson relation without altering the micro- 
scopic spectral correlators. This provides an indication that that lattice gauge 
simmulations can satisfy chiral random matrix statistics for finite lattice spac- 
ings even when genuine chiral symmetry has been sacrificed in the intrest of 
eliminating fermion doubling. 

Acknowledgements: The author gratefully acknowledge discussions with 
P. H. Damgaard, A. Andersen, J. Christiansen and A. D. Jackson. 
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